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ABSTRACT 


The results of coherence measurements of three free vrecession 
proton magnetometers are discussed. The measuring and data proces= 
sing instrumentation is described. The noise associated with the 
instrumentation has been measured and presented in statistical form. 
The results of measurements with elementary local perturbations are 
included to illustrate the target detection problem. The limitations 
of the instrument for target detection are discussed anid imsorvements 


suggested. 


The writer wishes to express appreciation to Professor Carl F. 
Menneken for suggesting the problem, and Professor Harold A. Titus 


for encouragement and assistance in the preparation of this thesis. 
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le Introduction 

The use of magnetometers for target detection requires identi- 
fication of signals in the presence of the noise associated with the 
earth's magnetic field. The signal, in many cases, will be smaller 
in magnitude than the normal diurnal variation of the field, Local 
fixed perturbations and geomagnetic pulsations. The system noise 
is of major concern in a detection scheme using magnetometers. This 
paper presents the preliminary results of a study of three proton 
free precession magnetometers operated under non laboratory condi- 
tions at Naval Postgraduate School La Mesa Village site, 

The free precession magnetometer consists of a coil of wire 
in which a source of protons is placed. A strong field is applied 
to orientate the proton spins in a direction approximately normal 
to the earth's field. The polarizing field is then removed ina man- 
ner such that the spins cannot follow the field. The earth's field 
then causes the protons to precess about the lines of magnetic force 
at the Lamour frequency given by: 

(UW = Xp H 

where Vp 2 gyro magnetic ratio 


magnetic moment/angular momentum (1) 


i 


and 


H magnitude of external field 

(earth's field) 

The net precessing magnetic moment causes a voltage to be induced in 
a coil surrounding the sample. Since the precessional frequency is 


dependent only on the external magnetic field, the measurement of H 


consists of measuring the frequency of the induced emf. The inherent 





accuracy of the free precession magnetometer is determined by the 
spectral line width (1), which is related to the relaxation time 
To, and is given bys 

| 


a a TC 


\ # 
Ye Tp 
where 3.6 ¥ (2) 
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vb 


H # Half intensity width 


T> 2 1 sec (meas) 


00 
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The gyromagnetic ratio has been measured by the NBS to tetter than 

1 part in 10°. Full advantage of this inherent accuracy has not been 
utilized because of the problems associated with the frequency measure= 
ment. The frequency measurement is complicated by the fact that only 
a limited time is available for measurement. The signal received 
from the magnetometer is in the few microvolt ranges hence locally 
generated electronic noise and transients from the polarizing opera] 
tion cause the signal to noise ratio to be small. Standard counting 
techniques are used for measurement. To get the required accuracy, 
the frequency must be obtained from period or multiple pericd average 
measurements. The process is digital in nature and lends itself 
very nicely to computer processing of the data. 

The main disadvantage is that any reasonable attempt to obtain 
an analog signal proportional to the measurement is an analog of 
period rather than frequency; therefore, the analog is the inverse 
of the quantity of interest. This does not present serious problems 
for target detection since the absolute magnitude of field is not 


of interest. 





2. Instrumentation 
2.1 Measurement 

The engineering aspect of the measurement problem can be stated 
simply as follows: 

Given: Three transient type signals of frequency approximately two 
ke/s occuring at the same time, with usable duration of about 3/) 
sece The signals are in the microvolt range. The recurrance rate 
is to be 2 sec. The measurements must be conducted in the presense 
of high electronic noise such as relay transients, contact noise 
and transients from the collapsing polarizing field. 

Find: The frequency of the oscillations accurate to 1 part in 10° 
or better, and suitable means of data aquisition and processing. 

The funds available excluded the possibility of using commer- 
cially constructed instruments or components. Also, modifications 
and interfacing would have to be done if components were used. Hence, 
the basic measuring instrumentation was constructed at NPGS. No 
detailed circuit design is presented; however the circuit diagrams 
are included in Appendix IV. 

Assume the precession signal has been amplified in a suitable, 
low noise, narrow band preamplifier, and converted to a signal have- 
ing fast waveform. Referring to the simplified block diagram of 
Fige 1, the basic counter operation is as follows: Let the slow 
and fast counters consist of flip-flop counting units, connected to 
form a binary counter. The gate control flip-flop is capable of 
being triggered into either state. The multiplier converts the clock 


frequency to a harmonic of the precision 100 kce/s oscillator. 
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Figure 1 

Let all counters be reset and ready for operation. The unknown 
input signal is applied to the gate. After a delay of approximately 
150 ms, an external signal triggers the gate control to open the sige 
nal gate. The long delay is necessary to allow for damping the polarie- 
zing field and relay transients. The first pulse through the signal 
gate opens the clock gate. The fast counter then counts the number 
of clock pulses occuring during the time interval determined by the 
slow counter. The slow counter is increased one for each cycle of 
unknown input signal. After the slow counter has reached a predetere 
mined number, both signal and clock gates are closed, and the number 
in the fast counter is a value of the time interval of 28 cycles of 
unknown input signal, or can be reduced to the average period of the 


signal. 





Let us define 
N 2 number in fast counter 
an 


Y 


t, 2 period of unknown signal 


60 


number of cycles of unknown signal 


period of clock ( W seco) 


Then we have 


TIn= Nt, 
t= lim 
IN a 

= Nop /2" 


Substituting in (1) gives 


. ae 
H Xp NTp 


and with n = 10, 


LY = 2348400 x /0o24 





(3) 


(4) 


The accuracy of such a scheme depends ons 

a The accuracy in determining the zero crossing of the unknown 
signal. 

be Propagation delay in the slow counter. 

Ce. Gating speed. 

d. Accuracy of the clock. 

e. Plus or mimus the las bit, due to noncoherent gating as 


illustrated in Fig. 2. 
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Figure 2 


The errors due to determining the zero crossing are usually the domi-e 
nant error in using this technique; however when averaged over 102) 
cycles they become small and the largest error is the plus or minus 
last bit, as shown. The accuracy in determining the zero crossing 


in the presence of noise (2) iss 


Vv 
t, se 
where ¢ 
t, 2 maximum error in trigger time in sec. 
Vn = peak to peak noise in volts. 
S e Slope of ideal signal at point of trigger level, 


For a sinusoidal inputs 


S = (Ns Sin we) 
= Vs W 
Assume a signal to noise ratio 2 V</Vp 2 20 dby and 
WE Zrxe Kos 
then “4 
t= Ons % 7/0 


w= IW TP 








For 102 pericdss 
tn® QO. 5s xe tag 
= 0.035 “4 5e 
The plus or minus last bit error if the clock is 100 ke, is +10 
S@Ce, and the other errors mentioned are negligible compared to the 


10 ) AY SCE 6 SILO» 


Consider the #1 bit limitation in detail for the earth's field 
50950 Oo » Which is the approximate value at Monterey, California. 
lets 

N = OOLZh051e e 1 bit 

1 bit #10 ) AY SEE « 
thens 


A Hel.) 4 


This becomes tie sensitivity of the instrument. If the clock - 


ed ° 
KA? 


increased to 1 mc/sec, 1 bit = Ae ikies ands 

A Hee Oil y gs however other noise becomes appreciable anid 
cannot be neglected, The circuits were designad for ] m/sec opera= 
tion, but the measurements to date were conducted at 100 ke/s and 
300 ke/s clock rates. The added complication is to construct 3 
meters that will read identical, the input circuitry and slow count- 
ers must have identical characteristics. This is impcssible fron 
a practical point of viex. The result is that the three counters 
with an ideal signal input can vary #1 bit in comparirg their rsad~ 
ings. 

A simplified block diagram and photo of the magnetometer is 


shown in Figs. 3, 4 and 5. It must be remembered thai operation 





is to be in remote locations where extremes of temperature, line 
voltage fluctuations, poor grounds, etc. will be founds; hence cir= 
cuits that are sensitive to these variations must be carefully come 
pensated. 

For the experiments under consideration, locating the sensing 
coils approximately 500' apart would be sufficients however for fur- 
ther experiments a distance of greater than 500! might be necessary. 
Therefore, a provision to place the preamps in the line nearer the 
sensing coils for improved signal to noise ratio was included. Bias 
for the preamps is applied to the signal line, from the console to 
eliminate the need for battery replacement. The relays that provide 
the two to three amps of polarizing current, and provide for dissapa-= 
tion of the stored energy when polarizing current is turned off, must 
precede the preamp. Since this unit is placed near the sense coil, 
it must be magnetically clean. 

Referring to the block diagram the components that are common 
to the three units are the clock, timer, recorders and the power 


supplies. 
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2.2 Data Aquisition and Processing 
Figure 6 shows a block diagram of the data collection and reproe 


auction system. 

The scheme used for data handling is certainly not optimun in 
the sense of tape utilization and ease of operation, but an Ampex 
CP1LOO was available and required no acditional expenditure fer tape 
recorder. Also with minor modifications a data recorder can be utie 
lized at a later date if funds besome available. 

Consider the form of the data. The sensing coils are polarized 
for one second, at the ani of this time the polarizing current is 
turned off, the transients are damped3 and the signal is fed to the 
counting units. After 102) cycles of the input signal (about 1/2 
sec.) the measurement is complete. In the fast counters, there are 
three binary numbers to be processed. The magnitude of the binary 
numbers, for H 2 50950 ¥ with clock frequency s 100 ke/s, is about 
OOL3W146, or 16 bits. With the clock increased to 300 ke/s, the 
number is 00h25260, or 18 bits. Hence, at the end of the counting 
interval, there are 48 bits to te recorded in a format for computer 
input. The Ampex CP100 is a "portable" instrumentation recorder 
with 1) parallel inputs and capable of speeds from 1 7/8 ips to 60 
ips. All of the record and reproduce anplifiers that were on hand 
were not capable of recording D.C. levels; therefore, tne information 
was put on tap2 in the form of vulses. Also, for reasonable speed 
in reproduction, the data was recorded at 1 7/8 ips and reproduced 
at 60 ips. A record of approximately 6.5 hours can be reproduced in 


10 minutes. The machine used to read the raw data is a CDC160 con=- 
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Figure 6 


cas, 





puter capable of reading words of up to 12 parallel bits. The mag= 
netometer words are broken into two, 12 bit words and time muitiplexed 
on tape to satisfy the computer input conditions. Hence, for each 
sample, there are six 12 bit words to ve processed. Even with the 
time speed up in reproduction the CDC 160 is operating below its 
capabilities. In an effort to utilize the 160 time, the data was 
put on tape, so that when it occured, the 160 would read at its maxi- 
mum rate and use the "no data time! to continue with preliminary 
processing. Channels 13 and 1 on the tape are input ready and dise- 
connect signals used to control the 160. And of course, for reason= 
able machine time efficiency, there is on hand enough tape to allow 
continuous recording for one week. 

The Ampex CP1CO was not designed specifically for this type of 
operation; hence the need for an interface to mate it with the 160, 
The interface must provide the following functions: 

a. Convert the pulse type data to D.C. level signals. 

be Put the data on line when a request is supplied from the 
computer. 

ce Effectively disconnect itself from the lines when other 
equipment is called. 

d. Provide eentnen information, i.e. input ready and discon- 
nect for the 160. 

e. Effectively isolate the computer from spurious transients 
from the tape machine. These normally occur during start <= stop 
operations, but the recorder is somewhat susceptible to powex line 


transients and they sometimes appear in its output. 


1h 





Figures 8, 9 and 10 shows the block diagram of the interface 
provided. Fig. 11 is a photo of the preliminary processing setup. 
The cards were assembled and mounted as a component on a laboratory 
cart, the tape unit is placed on the top of the cart and computer 
cable connectors provided. This eliminates the need for a permanent 
rack in the computer laboratory. 

The preliminary processing consists of the following: 

a. Reading the raw data. 

be. Checking the data and keeping a record of the errors. 

c. Writing the storage tape. 

d. Providing identifying blocks on the storage tape. 

e. Searching the storage tape for any particular record or 
the end of the data so that one storage tape may be used for many 
week's records. 

f£. Direct display of the raw data. 


The data storage format used is illustrated in Fig. 7. 


fe] ie 


STORAGE TAPE FORMAT 


22 BLOCKS 22 BLOCKS 
DATA DATA 


Figure 7 





With the time base used there are nearly 3800 twelve bit words 
per 6.5 hours record. The 6.5 hours record is one 3600 roll of ine 
strumentation tape at 1 7/8 ips. Since the operation is stopped to 


change tape, it is convenient to put the identifying block at the 
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end of this record. Also, this number of data words can be eracoae 
sed in the order of minutes by the CDC160l so that if the computer 

is needed for other jobss3 it can be stopped at the identifying blocks 
and return later to the same point with the proper index setting. 

The block length was set somewhat arbitrarily. The memory storage 
limit of the 160, write time of the 163 tape unit and suitable re- 
cords for handling in the 160), are the factors to be considered when 
chosing the block length. The instrumentation recorder cannot cone 
veniently be stopped when writing the storage tapes; hence one data 
point is lost when transferring data. This represents 22 pecints 

per 6.5 hours record or less than a minute in real time. The block 
length was chosen to be 356,9 twelve bit words per block, the fig- 
ure being compatable with the above conditions, 

The identifying record consists of writing the month, day and 
year in flex code. Tne remaining cells of the block are Loaded with 
zeros. It is then easy to separate the data from the identifying 
blocks since none of the data words will be all cgeros, 

The fact that the data had to be reconstructed in the 160), was 
also considered in chosing the block length. The information is 
buffered into the 160): in 6 bit words. The blocks consist of an 
even number of data points, i.e. six 12 bit or two 48 bit words. 
After a block is read into the machine and loaded into the A register, 
the magnetometer readings can be stored sequentially in memory with 
a loop of a few machine language instructions. The data analysis 
is then done using Fortran, allowing easy and versatile manipulation 
of the data. The 160 and 160); input programs are included in Appene 
dices II and III. 
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Reassembly is illustrated below. 
160, WORDS 


a 


Mag 1 (2h bits) Mag 2 (2h bits) 


Mag 3 (24 bits) 
Weg 9 (2h bite) | 


Briefly, the operation is as follows: Magnetometer readings 














are recorded for long periods of time at a remote site on a protable 
magnetic tape unit. The tape unit is then returned to the process- 
ing site, and using a time speed up of X 32, the intermediate proe- 
cessing is performed and the data stored. 

The storage tape along with suitable programs is then placed 
on file for processing by the computer facility at their convenience. 
The tape unit is then returned to magnetometer site for more records. 
The new data obtained is added to the storage tape so as not to dee 


stroy the existing data. 
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3. Data Analysis 

The experimental measurement program was set up to study the 
coherence between the magnetometers, i.e. the comparison of simulta- 
neous samples of the field, with the sensors located at different 
spatial locations. The coherence was to be examined with various 
local perturbations in the field. With the target detecticn objective 
in mind, the least complicated method of obtaining the target signal 
will result in the best instrumentation for field use. The problem 
is somewhat analogous to the radar problem where decisions must be 
made in the presence of noise, with associated probabilities of error 
in detection. 

There are several types of geomagnetic noisé that occur. The 
most dominant is the diurnal variations. Superimposed on this are 
smaller variations known as micropulsations and noise caused by lo= 
cal disturbances. Sudden eruptions on the sun sometimes cause large 
fluctuations in the field. Several studies of the various noises 
have been made (3), () and (5). For the proton precession magneto- 
meter, the instrument produces noise of major contribution to the 
spectra of small geomagnetic pulsations. For target detection, the 
signal of interest will most likely be of the same order of magni- 
tude as the small micropulsations and local distrubances,. For géoe= 
physical prospecting and mapping, the signals will be much iarger. 

As with any raw data, there will be some points that are obvie- 
ously in error. The cause of errors is usually uncertain, but some 
provision must be made in the processing to eliminate these errors 


So as not to distort the data. Care mst be taken when performing 
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these operations so the target signals are not masked. Records of 
the errors should be kept in order that decisions about the usefule 
ness of the set of data can be made. 

The problem of coherence between the spatial samples requires 
that the instrument noise be determined, with an ideal signal, then 
measurements be conducted with no local perturbations in the field, 
and the samples compared. The simplest comparison is to plot the 
difference signal vs. time. With no local perturbation, this means 
essentially a plot of the system noise vs. time3 hence the need for 
filtering is obvious. To use this method for target detection will 
require serious consideration of the filter time constants. 

The types of filtering 

1. <A clipping filter (on the raw data) that clipped at 4 200 Y. 
The number of points outside of this range are recorded. This range 
was about twice the expected diurnal variation. The filter removes 
the errors from lost high order bits and other instrument sources. 
The typical error rate was about 1/50000. 

2e A simple three point average for display of the raw data on 
the DD 65 display unit was used. This gives an operator a fast vise 
ual display of the raw data. Decisions can then be made about the 
quality of the data. 

3. Eighteen point averaging, i.e. 36 seconds real time, gave 
partially smoothed 6.5 hours graph records, 

he A simple single section low pass filter was used of the 
form 1/(s = a). The time constant was set at 100 sec. and 20 min. 


in real time. The implementation was done using Z transform (6). 


es 








The filter design was done by Professor Harold A. Titus. 


We have 


R(S) sg 


Xin(S) -- 





which has a Z transform 


x, 
X,.(2) Pinus eae 


In block diagram form this is scen in Fig. 29 





Figure 29 


In terms of the discrete time samples we have 


Xo(kti)= B = oe" 4, Lae 


M=O 
3 ge Ve eon 
neo 
T = sample period 
@ = time constant 


The function was trunkated at ns 20 for values of k S 20. 


The contribution to the filtered value from the terms multiplied 


by e239 for n ‘> 20 is negligible. 


2h 





The coherency is a statistical parameter, and must be repre- 


sented as suche 


Arth- %® 
A\ 2 
A 3# %2 > %3 


The differences of the unfiltered values were grouped in a frequency 


Let: 


Af 


I = H3 


distribution which approximates the continuous density function for 
large sample size. Since the process is digital, the allowable dif- 
ferences, in bits ares 0, 41, 4 2, $ ceo, Therefore, the frequency 
distribution is not a smooth curwe.e Since the instrument noise is 
randomly distributed, the envelope of the distribution is a normal 
density function. The mean was removed before plotting, and the 


statistic representing the standard deviation was computed according 


tos 
== 2. 
~ | F(Ai-A 
4! 
where s N 2 the number of samples in a 6.5 hour record. 


Comparison of the density function curve with an ideal signal 
input to that of magnetometer input gives a measure of the noise 


produced by the device. 


To further define the coherence, the auto and cross correla= 


tion functions were computed and the coherency computed as defined: 
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This discrete computational form is 
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N 2 number of data values 
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nunber of lags 
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Since the functions are from finite time series n,,, can only be made 


about 0.1 N. The 6.5 hour records will vary in form deverding upon 


the time of day the records wers mace. Therefore, the correiation 





functions will not be the same shape for each record; however the 
coherence function as defined will be the same, and for perfectly 
coherent records is unity. The programs are included in Appendix 
Til. 

It should be pointed out that when perturbations are placed in 
the field, the basic assumptions regarding correlation; i.€. a sta= 
tionary random process, are no longer applicable. 

Experimental data was taken to determine typical signatures for 
elementary perturbations. Fig. 12 shows the geometry of the detector 


centered coordinate system. 
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The dipole moment was orientated parallel to the earth's field and 
its motion was in the x direction such that dx/dt e constant. The 


approximation r )>1 is valid. 


f= \ Ho” = Ho* 


At the detectors 


And3 
Hp 3 KM 
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dipole moment 


function of time and the geometry 
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No attempt was made at this time to compare the measured values with 
the theoretical values, because the terrain does not lend itself 
keeping L, d, and v constant. The results are representative of 


typical signatures that would be found in practice. 
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4. Results 

hel Instrumentation: 

An ideal signal was substituted for the procession signal and 
set of records were made, as shown in Figs. 13 through 25. The 
graph scales were left. unchanged because it was desired to show the 
instrument errors on the same scale as the magnetic field scales, 
The magnetic field vs. time graphs were not included because the in- 
put signal was an audio oscillator and only 4 small drift puts it 
off the graph scale. The clock was 300 ke/s which makes 1 bit s 
0.36 ¥ . 

Figure 28 shows the precesion signal from the various sensors. 
Figs. 26 and 27 are typical plots made from the analog section of the 
instrument. For comparison, the Rb87 magnetometer outout was adjust- 
ed to the same scale and plotted on the same graph. The nb 87 magnetoe 
meter output is an analog of frequency, while the proton magnetometer 
the analog is of period; hence the reason for the opposite drift on 
the chart. Figure 27 is an analog record of the differences. The 
filter was a single section RC filter with time constant = 10 sec. 

On the multiple curve plots showing coherence, the curves are identiec 
cal and cannot be distinguished from each other. These curves can 
be compared with the curves obtained from the sensor with no local 
perturbations in the field, for the contribution of the measuring 
instrument noise to the total noise. 

Difficulties with both instrumentation and computer tape have 
been encountered. The use of new computer tape for data storage and 
scratch tapes, and frequent head cleaning of the instrumentation 
recorder reduced the difficulties to a minimum. 
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Differences in magnetic intensity vs. time 
for an ideal signal input. 

DELL = Ky = Ho 

DEL2 = Hl - H3 

DEL3 = Ho = H 

Filter a 18 point average 
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Differences in magnetic intensity vs. 
time for an ideal signal input. 


DEL 2 = ki ~ H3 
DEL 2 regio = H3 
Filter = Z transform low pass with 
time constant of 20 minutes. 
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Normalized density function for the 
difference HH) - Ho - mean (H] - Ho), = 
with an ideal signal input. © 
Mean = 0.061 gamma 

Standard deviation = 0.21 
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Normalized density function for the 


difference Hy - H3 = mean (Hy - H3), = 
with an ideal signal input. © 
Mean 2 0.056 gamma 
Standard deviation = 0.72 
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Figure 16 
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Normalized density function for the 


difference Ho - H3 - mean (He - H3)s = 
with an ‘deat econ input 
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Standard deviation = 0.72 
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Auto correlation for an ide2l 
Signal input. 


t Record length = 6.5 hours 
8 Hy = Correlation (fy = meen 4j) 
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Cross correlation for an ideal 
Signal input. 





f Record length = 6.5 hours 
* Hig = Correlation (Hj - mean H) and Ho - mean Ho) 
H13 = Correlation (H. = mean Hy and H3 = mean H3) 
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Figure 19 
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Coherence for an ideal signal inout 
Record Length = 6.5 hours 

Hio = Coherence of Hy and Ho 

H13 = Coherence of Hy and H3 

Ho3 = Coherence of Ho and H3 
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Differences in magnetic intensity 
for an ideal signal input 

Record Length = 1 hour 

Filter = 2 point average 
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Difference in magnetic intensity 
for an ideal signal input 


§ Record length = 1 hour 
Filter = 2 transform low pass 
with time constant of 100 seconds. 
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Auto correlation for an ideal 
Signal input. 

Record length = 1 hour 

Hy = Correlation (H | - mean H)) 
Ho = Correlation (Ho = mean Ho) 
H3 = Correlation (Hp = mean H3) 
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Figure 23 
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Record length = 1 hour 
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Coherence for an ideal signal inout 
Record length = 1 hour 

H1o = Coherence of Hy and Ho 

H13 = Coherence of H, and H3 

H23 = Coherence of Ho and H3 
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TYPICAL ANALOG PLOT OF 
THE DIFFERENCES VS. TIME 
FILTER = RC WITH TIME CONSTANT OF 6 SECONDS 
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h.e2 Magnetic field measurements 





Only representative results of the magnetic field measuronents 
are presented, since the measurements are not complete at this time. 
Samples of plots made with the clock at 100 ke/s are included to 
illustrate the cffect of the clock frequency on the results. For 
the associated density function of the cifferences, the allowable 
values are 021.1, E 2,2, So © © in gamma, with the clock at 100 
kc/se The expanded difference vs. time shows clearly the ¥ last 
bit limitation. Comparing the density functions with the plots for 
the clock at 300 ke, shows that most of the noise associated with 
the measurement is due to the Last bit and can be reduced by increas- 
ing the clock frequency. The perturbation experiments were conducted 
with the clock at 300 kc/s. The resutts of the perturpation experie 
ments show clearly the effect of the digital filtering, and the ime 
portance of tame constant consideration for target detection. Thea 
oscillator is to be run 1 mc/s for further experiments. Also, many 
of the errors that appear in the data were the direct result of the 
instrumentation tape. New tape has beén purchased and should elim= 
inate most of these errors, in future records. 

As pointed out oreviously, the basic assumptions regarding core 
relation are not valid, when the noise is not stationary; however, 
since correlation detectors can be instrumented and may vrove useful 
in the detection problem. Some of the correlation function plots are 
included to illustrate the effect of local perturbations on the funce 


tions, 





S. Conclusions 

For target detection, where the signals under consideration are 
in the 4 gamma, or less range, it will be difficult to detect the 
target without resorting to sophisticated signal processing. In 
this range of signals, consideration should be given to other instrue 
ments; such as the Rubidium ootically pumped magnetometer, whose 
theoretical sensitivity is 0.O0L gamma, especially where the detector 
is stationary, and orientation is nota problem. It can be demonstrated 
that wien using the same measurement technique as for the proton 
magnetometer, the instrument and processing noise will be an order of 
magnitude less than the theoretical limits imposed by the sensiti-e 
vity of the instrument. 

For signals in the $ to 14 gamma range, detection could be ace 
complished with analoz type correlation detectors. These could be 
implemented with operational amplifiers. 

Signals above this ran je could be detected with the simple dife 
ferencing technique, and the associated R.C. filters. 

Increasing the clock frequency and averaginz the period over 
more than 102) cycles would reduce the measurement noise to a smaller 
part of the total noise. The target signals could then be reduced 
by a factor of three, with no change in the instrumentation. 

Full utilization of the inherently narrow line width of proton 
magnetometer scannct de achieved using the starlard measuring tech 
nique. Other techniques have been suggested by Professcr Carl I. 
Mennekens such as utilizing a rapid soherent polarization in a phase-= 
lock loop arrangement. Further devolopment work needs to be done 


to explore these possibilities. 
hy? 
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Difference in magnetic intensity vs. 
time for clock at 300 kc/s and no 
local perturbations. 

Record length - 5.5 hours 

DEL 1 = 4, - Ho 

Dsl 2 Ry = a 

DEL 3 = Ho - 43 

Filter = 18 point average 
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Pa Difference in magnetic intensity vs. 
& time for clock at 300 kc/s and no 
local perturbations. 

Record length = 6.5 hours 

DEL 2 = - H3 

DEL 3 eee > H3 

Filter - Z transform low pass with 
time change = 20 minutes 
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Normalized density function of the 
difference,D. No local perturbations 


Record length = 6.5 hours mH 
Clock = 300 ke/s S 
(] ~ Ha) - mean (i, = Hp) =D 
Mean = -.09 gamma 
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Normalized density function of the 
difference,D. lb local perturbations 
Record Length = 6.5 hours 

Clock = 300 ke/s 

(Ho - H3) - a oe = op =D 
Mean = 25, 33 ¢g : : 
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Normalized density function of the 
difference,D. No local verturbations 
Record eee - 6.5 hours 

Clock = 300 ke/s 

(Hy - H3) = mean (Hy - 43) =D 
Mean = 3. 70 gamma 
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5 Auto correlaticn for earth's magnetic 

field vs. time. No local perturbations 

Record length = 6.5 hours 

Clock = 300 ke/s 

Hy = Correlation (Hy, = mean dj) 
Wy Ho = Correlation (Ho = mean Ho) 
& H3 2 Correlation (H3 - mean H3) 

lag = At 
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> field vs. time. No local perturbations 
Record length = 6.5 hours 
Clock = 300 kc/s 
Ho = Correlation (4) - mean Hy and Hp - mean Ho) 
Hi3 = Correlation (Hj = mean Hy and H3 - mean H3) 
\g Ho3 - Correlation (Ho - mean Hy and H3 = mean H3) 
8 Lag = At 
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Coherence for earth's magnetic field 
vs. time. No local perturbations 
Record length = 6.5 hours 
Clock = 300 kc/s 

_Hyo = Coherence of H, and Ho 
H13 = Coherence of th and H3 
H23 = Coherence of Ho and H3 
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Barth's magnetic field vs. time with 
2O gamma step function applied to H3¢ 
Record leneth - 6.5 hours 

Clock = 300 ke/s 

Filter - 15 point average 

Mean = (mean Hy + mean Ho + mean H3)/3 
Dies, hye = eat 

io = Hy = Mean 

13 = H3 - Mean 

Mean = 0.509307 E + 05 
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Barth's magnetic field vs. time with 
20 gamma step function anplied to He 
Record length = 6.5 hours 

Clock = 300 ke/s 

Filter = 2 transform low pass with 


053 


time constant of 20 minutes 

Mean = (mean Hy) + mean Ho 4 mean H3)/3 
Ay = By - Mean 

Ho = Ho - Mean 

H3 2 Hq - Mean 

Mean ~~ 0.509307 E + 05 
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Difference in magnetic intensity vs. time 
with 20 gamma step function applied to H3- 
Record length = 6.5 hours 
Clock = 200 ke/s 
Filter = 18 point average 
DEL L s Hy ~ Ho 
DEL 2s hj) - H 
DEL 3 - H» = i 
nd 
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Difference in magnetic intensity vs. time 
with 20 gamma step function applied to H3. 
Am Record length = 6.5 hours 
S Clock = 300 ke/s 
Filter = Z transform low pass with 
time constant of 20 minutes 
DEL 1 = - Ho 
DEL ¢ = Hy = H3 
DEL 3 = He - H 
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Barth's magnetic field vs. time 
Filter = 18 point average 

Clock = 300 ke/s 

Record Length s 6.5 hours 

Hy = H) = mean hy 

Ho = Ho - mean Hy 

Ha 2 He - mean Fy 


Perturbations: 

+ 5 gamma and + 10 gamma applied 

at constant velocity to H2. = 2 
gamma, - 5 gamma, - 5 gamma applied 
at constant velocity to Hn, Ho and 
H, curing first hour. A 5 gamna 
step applied to Hy at t = 1 hour and 
removed at t = 1 3/l hours. 
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2 Barth's magnetic field vs. time 
oS Filter = Z transform low pass with 
time constant of 20 minutes 
Clock = 300 ke/s 
Record length - 6.5 hours 
Byes Hy = mean } 
, oo - yt - mean 2 
S 13 = Hy - mean Hy 
Perturbations; 
+ 5 gaia and + 19 gamma applied 
at constant velocity to Ha. - 2 
gamma, ~ 5 gamma, - 5 gamma applied 
12 at constant velocity to H,, Ho and 
H, during first pones A > garma 
step applied to H t - 1 hour and 


removed att -1 3h, hours. 
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Earth's magnetic field vs. time 
Filter = 2 point average 

Clock - 300 ke/s 

Record length = 1 hour 


Hy = Hy - mean Hy 
Ho ~ Hy - mean Hy 
H2 = Oo, - méan Hy 


Perturbations: 

+ 5 gamma and + 10 gamma applied 

at constant velocity to H.. 2 
gamma, -5 gamma, -5 gamma” applied at 
constant velocity to Hz, Ho and Hy 
during first hour. A 5 gamma step 
applied to iy at t = 1 hour and 
removed at t= 1 3/\; hours. 
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Earth's magnetic field vs. time 
Filter = 2 transforn low pass with 
time constant of 100 seconds. ° 
Clock = 300 ke/s 

Record length = 1 hour 

Hy = Hy ~ mean H 

Ho = io - mean Hy 

13 = H2 - mean Hy 

Perturbations: 

= 5 gamma and = 10 gamma applied 
at constant volocity to 44. - 2 
gamma, - 5 gamma, - 5 gamma applied 
avy constant velocity to H » H, and 
Hy during first hour. A 5 garima 
step applied to H4 at t = 1 hour 
and removed at t 21 3/) hours. 
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Hj 
Auto correlation for earth's magnetic 
merela VS. time. 
Record Leneth x 1 hour 
Clock = 300 kc/s 
Hy = Correlation (H, = mean H}) 
Ha = Correlation (Hj - mean H) 
H3 = Correlation (H3 - mean H3) 
Lag = At 
Perturbations: 
4 5 gamma and + 10 gamma applied at 
constant velocity to H.. -2 gamma, 
-5 gamma, -5 gamma applied at constant 
velocity to N35 Ho and Hy during first 
hour. A 5 gamma step applied to He at 
% = 1 hour and removed at t = 1 3/h hours. 
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HOI 
H2 1 

Cross correlation for earth's magnetic 
field vs. time. Hi3 


Record length - l hour 

Clock = 200 ke/s 

Hy5 = Correlation (H, -mean H, and H, -mean Ho) 
H13 = Correlation (Hy -mean Hy, and H3z -mean H3) 
Ho3 = Correlation (Hp -mean H», and H3 -mean H3) 
Bags A © 

Perturbations: 

+ 5 gamma and 4 10 gamma applied at 

constant velocity to Hy. -2 gamma, 

-5 gamma, -5 gamma applied at constant 
velocity to Hz, Ho and H, during first 

hour. A S gamma Step applied to H, at 

4 = L hour and removed at t = 1 3/f hours. 
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Coherence for earth's magnetic field 
vs. time. 

Record length - 1 hour 

Clock = 300 ke/s 

Hy5 = Coherence of Hy and Ho 

Hy3 = Coherence of Hy and H 

a3 = Coherence of Ho and H 

Lag - At 

Perturbations; 

+ 5 gamma and +10 gamma applied at 
constant velocity to 4... -2 gamma, 

-5 gamma, -5 gamma applied at constant 
velocity to H3, H, and Hy, during first 
hour. A 5 gamma step applied to H, at 
t = Ll hour and removed at t - 1 3/t, hours. 
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160 PROGRAMS AND OPERATOR INSTRUCTIONS 


Starting Addresses 


‘Q000 
0002 
0004 
0006 
0010 
0100 


Error Stops 
ERRO1 
ERRO2 


ERRO4 
ERRO 3 


Addresses of 
Important Counters 


0057 
0045 
0046 


Function 
Rewind 
Read 
Search 
Write 


Write identification Block 
Main Input Frogram 


Cause 


17g of current input in error 
Memory storage net complete 
End of storage tape 

163 not ready 


Information 


Number of blocks written on 
storage tape 

Current storage location in 
memory 

Number of zero words in meme 
ory 


Operator Instructions 


1. Load program and make the necessary cable connections, 


2. Set identifying records in cells 30g to 47g,. 


3. Write identifying block (run from P s 010). 


ir = 100. 
4, Set cell 100 = 


= 


Si s 


5. Start Ampex C.P. 100 and run from P z= 100. 


2200 
O 


stop at 


Stop when 


tape is finished or ERR stop. 


II-=-1 





° 
S 
= 


ie 7 th 
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INPUT PROGRAM 


ERRO! 


Tage 
Ci 








AOD 
CNTR 
NO 





Il-2 





0100 
0101 
0102 
0103 
0104 
0105 
0106 
0107 
0110 
0111 
0112 
0113 
0114 
0115 
0116 
0117 
0120 
0121 
0122 
0123 
0124 
0125 
0126 
0127 
0130 
0131 
0132 
0133 
0134 
O15 
0136 
0137 
0140 
0141 
0142 
0143 
0144 
0145 
0146 
0147 
0150 
0151 
er52 
O55 
0154 
O15 
0156 
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2200 
1000 
445 
0400 
4046 
0400 
4150 
4151 
4152 
4153 
A154 
4155 
7500 
0500 
(203 
0620 
6102 
0600 
1600 
0607 
6115 
2150 
6013 
2151 
6011 
2152 
6007 
2153 
6005 
2154 
6003 
2155 
6105 
5446 
0777 
6536 
0001 
2050 
4047 
2147 
ALAS 
5445 
5447 
3456 
6505 
2200 
ATTT 


IIl=3 


LDFOO 

1000 
STD45 
LDNOO 
STD46 
LDNOO 
STI50 
STI5i 
STI52 
STI53 
STI54 
STI55 
EXFOO 
0500 
INPO? 
0620 
NZFO2 

0600 
LSFOO 
ADNO7 
NZF15 
LDI50 
ZIF13 
LDI51 
yeas Vaal 
LDI52 
ZIFOT 
LDI53 
ZIFOS 
LDI54 
ZIFO% 
LDI55 
NZFO5 
AODA6 
SBN77 
NZB36 
ERROI 
LBD5O 
ST D47 
LDI47 
STI45 
AOD4S 
AOD47 
SBD56 
NWZBO5 
LDFOO 

4777 


(ADDR ) 


(CNTR1) 





a: a -_—o a 


ORESY 
0160 
0161 
0162 
0163 
0164 
0165 
0166 
0167 
0170 
Cale 


0044 
0045 
0046 
0047 
0050 
0051 
0052 
0053 
0054 
0055 


0057 
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3445 
6653 
2044 
BsaD 
6003 
6657 
0002 
7101 
0172 
1000 
7700 


ADDITIONAL CELLS USED 


LIMIT 

ADDR 

CNTR 1 
CNTR 2 
BUFFER 
BUFFER 
BUFFER 
BUFFER 
BUFFER 
BUFFER 


NUMBER 


STORAGE 
STORAGE 
STORAGE 
STORAGE 
STORAGE 
STORAGE 


SBD45 
PJB53 
LDD44 
SBD45 
ZIFO3 
PJB57 
ERRO2 
JFIOL 
WRITE 


HLTOO 


(PUT IN LIMIT) 
7600 


OF MEMORY DUMPS 
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MEMORY DUMP TO M.T. 
Se 







Ati’ 


Sok 
1& 
Appr 


Be 


Pa 
sf 

~J 

™~ 


Yes 


STORE 


© 
=i 


Str 
1ADDR 
SMD 
N 
AOD 
ADDR 
ie 1} ADDR 
us ) [taip 


a 


@ 


YES 


PUNCH 


II-5 








0172 
0173 
0174 
Ol(> 
0176 
0177 
0200 
0201 
6202 
0202 
0204 
0205 
0206 
O20, 
0210 
0211 
O2le 
S215 
0214 
Sea 
0216 
Oal7 
0220 
Gan) 
222 
229 


0163 
0164 
0165 
0166 
0167 
0170 
Cnty 1 
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Dato 
2600 
TTTE 
4060 
7500 
11421 
7600 
0202 
6004 
5460 
6506 
0003 
%500 
eae 
7 303 
7600 
6102 
LOOO 
7500 
1141 
7600 
0240 
6103 
(108 
O100 
7700 


DUMP TO PAPER TAPE 


{500 
4104 
(909 
7600 
6102 
1000 
1100 


TERM 


AQOD57 
LGFOO 


STD60 
EXFOO 


REQ STAT 


INAQO 
LPNO2 
ZIFO4 
AOQDEO 
NZBOG 
ERRO 3 
HAFCO 
ene 
OUTO 3 
OUT 
NZFO2 


Cor ieee 


TET WRITE 


EXFOO 


REQ STAT 


INAOQO 
LPN40 
NZFO 3 
JFIOL 

OO 
HLTOO 
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DAILY IDENTIFICATION 


II=-7 


USE 8 CELLS EXAMPLE 
0030 0070 & day 
0031 0064 9 
0032 0032 J 
0033 0030 A month 
0034 0060 N 
0035 0072 6 
0036 0062 Besse 
0037 0045 CAR RET 
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DAILY RECORD TO M.T. 
aoe VEU SU MOT 














0224 
0225 
0226 
S227 
0230 
Ozer 
0232 
0253 
0234 
0245 
0236 
0257 
0240 
0241 
0242 
0243 
0244 
0245 
0246 
O247 
0250 
O2bs 
0252 
0255 
0254 
0255 
0256 
0257 
0260 
0261 
0262 
0263 
0264 
0265 
0266 
0267 
9270 
0271 
Deve 
O27 > 
0274 
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DAILY RECORD TO M.T. 


2030 
4100 
LOOO 
2Oou 
4100 
LOOL1 
2032 
4100 
L002 
2033 
4100 
L003 
2034 
4100 
1004 
2035 
4100 
1005 
2036 
4100 
1006 
2037 
4100 
1007 
2200 
(700 
4100 
O100 
2200 
1010 
O45 
0400 
4145 
2044 
3445 
6003 
5445 
6506 
Leu 
On f2 
£100 


Lie 


LDD3O 
oTIOO 
LQOQO 
LDD31 
STIOO 
LOOL 
LDD3e2 
STIOO 
L002 
LpD>> 
ST LOO 
1003 
LDD34 
STIOO 
1004 
LDD35 
oT Ioo 
OOD 
LDD36 
ST IOO 
1006 
LDD57 
ST IOO 
LOOT 
LDFOO 
~700 
STIOO 
OLO00 
LDFOO 
1010 
STD45 
LDNOO 
SrI45 
LDD44 
SBD45 
2JFO3 
AOD4S 
NZBO6 
JFIOL 
Onl Ge 


(ADDR } 


(ADDR ) 


(DUMP ) 


NORMAL STOP 






NO 


REWIND, READ AND 
DAILY RECORD SEARCH 


LDN 
OO | 





REwIind 
ENARM 





er 





MALL s&h 
Ree EE oe 


REWIND, READ AND 


DAILY RECORD SEARCH 






START 
SEARCH 


ALPHA 


YES 


e a -_— 2 








Start 


Rewind 0275 


Read 
Start 


Search 


Start 


0276 
0277 
0300 
0301 
0302 
0303 
0304 
0305 
0306 
O307 
O3510 
O311 
O3l2 
0313 
0314 
CoU5 
0316 
O317 
0320 
0321 
O322 
0323 
0324 
0525 
0326 
0327 
0330 
O)oyoul 
O52 
0333 
0334 
0335 
0336 
O557T 
0340 
0341 
0342 
0343 
0344 
0345 
0346 
O347 
0350 
0351 
0352 
O35) 
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0400 
4043 
7500 
1141 
7600 
0202 
6504 
7500 
aod: 
7500 
1141 
7600 
0202 
6504 
7700 
7500 
1141 
7600 
0202 
6504 
T1500 
2111 
7203 
7600 
6102 
1000 
7500 
1141 
7600 
0240 
6002 
0004 
2043 
0701 
6004 
6525 
O401 
4043 
2100 
1000 
3430 
6532 
2100 
1001 
3431 
6536 
2100 


II-le 


LDNOO 
STD43 
EXFOO 
REQ STAT 
INAOO 
LPNO2 
NZBO4 
EXFOO 
REWIND 
EXFOO 
REQ STAT 
INAOO 
LPNO2 
NZBO4 
HLTOO 
EXFOO 
REQ STAT 
INAOC 
LPNO2 
NZBO4 
EXFOO 
EXF11 READ 
OUTO% 
TERM READ 
NZFO2 
INIT READ 
EXFOO 
REQ STAT 
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